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This paper presents the results of statistical analyses performed to predict the thrust imbalance between 
two solid rocket motor boosters to be used on the Space Launch System (SLS) vehicle. Two legacy internal 
ballistics codes developed for the Space Shuttle program were coupled with a Monte Carlo analysis code to 
determine a thrust imbalance envelope for the SLS vehicle based on the performance of 1000 motor pairs. 
Thirty three variables which could impact the performance of the motors during the ignition transient and 
thirty eight variables which could impact the performance of the motors during steady state operation of the 
motor were identified and treated as statistical variables for the analyses. The effects of motor to motor 
variation as well as variations between motors of a single pair were included in the analyses. The statistical 
variations of the variables were defined based on data provided by NASA’s Marshall Space Flight Center for 
the upgraded five segment booster and from the Space Shuttle booster when appropriate. The results 
obtained for the statistical envelope are compared with the design specification thrust imbalance limits for 
the SLS launch vehicle. 


I. Introduction 

T he knowledge of the difference in thrust at any point in time between two motors of pair operating on 
the same vehicle, referred to in this work as thrust imbalance, is necessary for predicting the loads on 
the vehicle and launch pad during motor ignition and for vehicle loads and control requirements during 
flight. This work uses a Monte Carlo statistical method approach developed and verified for the Space 
Shuttle boosters by Sforzini and Foster in Refs. 1 and 2 to develop statistical envelopes for the thrust 
imbalance versus time during both ignition transient and steady state motor operation of the Space 
Launch System (SLS) vehicle. The analysis is done in two parts using two different computer codes, one 
for the ignition transient and one for the steady state operation of the motors. The code chosen for the 
ignition transient was developed by Caveny and Kuo in Refs. 3 and 4 and the code used for the steady 
state operation was developed in Refs. 1 and 2. The codes were chosen because of their proven ability to 
simulate the internal ballistic performance characteristics in an accurate and efficient manner. The need 
to analyze a large number of motors, of the order of 2,000 individual motors in a reasonable amount of 
time made these two codes logical choices. The details of the internal ballistics codes used are well 
documented in the above cited references and are not being repeated in this paper. A MATLAB shell 
program was created and used to coordinate the user input files, program execution, output text files and 
plotting routines. The use of the programs in the MATLAB environment is described in Ref. 5. 
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II. Ignition Transient Thrust Imbalance Analysis 


The Monte Carlo method described in Ref. 1 was coupled with the one-dimensional internal ballistics ignition 
transient code developed in Refs. 3 and 4. This code was chosen since it is still currently being used by 
NASA/MSFC personnel and other organizations for ignition transient predictions and because it was readily 
adaptable to the Monte Carlo technique. A more complete discussion of the utilization of the computer code will be 
given later in this section. 

The Monte Carlo analysis of the ignition transient thrust imbalance used 33 selected variables for the ignition 
transient code. These variables are presented in in Table II. 1 below. 


Table II. 1 Monte Carlo Analysis Ignition Transient Variables. 


Variable Name 

Variable Definition 

TPI 

initial propellant grain temperature 

AT 

nozzle throat area 

XP 

location where propellant begins 

XG 

aft end of propellant 

XE 

end of flow passage 

GAMA 

ratio of specific heats 

W 

molecular weight of combustion gases 

TIGN 

igniter gas mean temperature 

TFREF 

reference adiabatic flame temperature 

RUFSUR 

port wall roughness 

FKPR 

propellant thermal conductivity 

ROPR 

propellant density 

CPR 

propellant specific heat 

TOREF 

reference propellant temperature 

SIGP 

temperature sensitivity of burning rate a constant pressure 

TPSCRI 

surface ignition temperature 

RREF 

reference burning rate 

PREF 

propellant property reference temperature 

BREXP 

burning rate exponent 

EBC 

Robillard-Lenoir constant 

EBEX 

Robillard-Lenoir constant 

DE 

diameter of the nozzle exit plane 

CM 

nozzle thrust loss coefficient 

EROAT 

nozzle erosion parameter 

EROEXP 

nozzle erosion parameter 

ALFAD 

nozzle divergence half-angle 

TPSHFT 

reference PMBT difference between motors of a single pair 

BRSHFT 

reference burning rate shift between motors of a single pair 

ATSHFT 

reference throat area shift between motors of a single pair 

DPSHFT' 

reference propellant density shift between motors of a single pair 

DESHFT 

reference exit diameter shift between motors of a single pair 

ETSHFT 

reference throat erosion rate shift between motors of a single pair 

SHFTMI 

reference At time shift in mass addition versus time from the igniter 
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Some of the variables listed in Table II. 1 were held constant for the present analysis, but can be used if their 
statistical variation properties are known. As can be seen from the variable definitions above, there are motor to 
motor variations and variations within a single pair of motors. Variations within a pair are identified by SHFT either 
before or after the variable name. For example BRSHFT is the within pair variation of the burn rate RREF. The 
analysis consisted of selecting input parameters for 2000 5-segment SLS solid rocket motors or 1000 motor pairs. 
The variations between motors of a pair were defined for what was referred to as the odd and even numbered motors 
in that pair. Specifically, the odd numbered motor was considered the base for the pair and the motor to motor 
variations were imposed on the even numbered motor relative to the odd numbered motor. Based on the Monte 
Carlo variable selection process, input files for the 2000 individual motors were created. These input files were then 
analyzed sequentially and individually in the internal ballistics ignition transient code in the internal ballistics 
ignition transient code in the internal ballistics ignition transient code to determine the thrust time trace for each 
motor. The results for the odd and even motors were stored separately and the thrust imbalance versus time was 
calculated for each pair of motors. 

The results presented below are for a propellant mean bulk temperature (PMBT) selected randomly between 
40°F and 90° F and a within pair variation given by TGSHFT. The burn rate reference temperature was 60° F. The 
thrust imbalance versus time data was calculated for each pair and the thrust imbalance for analyzed statistically at 
selected time points. The number of time points selected was the results of this analysis were the mean and standard 
deviation for the thrust imbalance at each time slice chosen. Based on this data a prediction of the K-sigma limit 
about a zero mean value envelope was constructed as a function time. For 1000 pairs, K was chosen to be 3.4725 5 . 
This value of K is for a probability of 99% that 99.9% will be included between the mean value and +/- K times the 
standard deviation. This envelope was then compared to the design thrust imbalance limits for a pair of motors 
during the ignition transient as shown in Fig. II.l. 
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Thrust Imbalance Envelope for 1000 Motor Pairs 



Figure II-l. Steady state thrust imbalance envelopes 
(PMBT mean = 40° - 90° F and K=3.472). 


The thrust imbalance data was also analyzed to determine the mean value and standard deviation of the 
maximum absolute value of the thrust imbalance and the time at which it occurred for specified time intervals during 
the ignition transient. In addition to the plot of the imbalance envelope, the results for all of the motors were 
available for plotting, both together and individually. Fig. II.2 presents a plot of plus or minus the absolute value of 
the thrust imbalance versus time for 1000 motor pairs compared to the envelopes shown in Fig. II.l. The results 
shown above indicate that for the statistical variations expected in the selected internal ballistic performance 
parameters that the resulting thrust imbalance between the two motors will meet the current design criteria limits. 
Furthermore, the thrust imbalance envelopes predicted by the Monte Carlo analysis agree quite well with the current 
design envelope with regard to the time dependent variation in thrust imbalance during the ignition transient for the 
SLS motors. It should be pointed out that the NASA thrust imbalance envelopes are valid for any odd motor of a 
pair having a PMBT between 40° F and 90° F with the even motor being constrained to be within the allowable 
temperature difference between motors of a single pair. We point this out because the NASA envelopes shown later 
in this paper for the steady state thrust imbalance have been scaled to eliminate the effects of PMBT on the thrust 
imbalance versus time. We will say more about this later in the paper. 
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+/- Absolute Value of Thrust Imbalance for 1000 Motor Pairs 



Figure II.2. +/- Absolute value of the thrust imbalance during the ignition transient 
(PMBT mean = 40° - 90° F and K=3.472). 

Fig. II.3 shows the thrust time traces for 2000 motors during the ignition transient. Figs. II.4 and II. 5 show the 
thrust time traces for the odd and even numbered motors respectively. These figures show the variability in thrust 
between the individual motors. To illustrate one feature of the possible data reduction capability contained in the 
codes Fig. II. 6 shows a typical thrust imbalance trace for a single motor pair compared to the thrust imbalance 
envelopes and Fig. II. 7 compares the thrust time traces for the same pair of motors. These latter two plots can be 
useful because they allows the analyst to locate motor pairs where the thrust imbalance may be high and then to look 
at the performance parameters for that particular pair to evaluate what may have caused the high imbalance. 
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Thrust Thrusl 


Ignition Transient Thrust for 1000 Pairs of Rocket Motors 



Time 

Figure II.3. Thrust time traces for 2000 motors during ignition transient 
(PMBT mean = 40° - 90° F) 

Ignition Transient Thrust for Odd Numbered Rocket Motors 



Figure II.4. Thrust time traces for odd numbered motors 
(PMBT mean = 40° - 90° F). 
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Ignition Transient Thrust for Even Numbered Rocket Motors 



Figure II.5. Thrust time traces for even numbered motors 
(PMBT mean = 40° - 90° F). 


Ignition Transient Thrust Imbalance for Motor Pair 123 



Figure II.6. +/- Absolute value of the thrust imbalance for a single motor pair 
(PMBT mean = 40° - 90° F and K=3.472). 
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Ignition Transient Thrust for Motor Pair 123 



Time 

Figure II.7. Thrust time traces for a single motor pair 
(PMBT mean = 40° - 90° F). 


III. Steady State and Tail off Thrust Imbalance Analysis 

The Monte Carlo method described in Ref. 1 was previously used and validated with regard to 
predicting the thrust imbalance for the Space Shuttle SRM’s, as discussed in Refs. 2 and 6 through 17. 

The Monte Carlo analysis of the RSRMV thrust imbalance during steady state operation used 
thirty six selected input variables for the ignition transient code. These variables are presented in 
Table IV. 1 below. 

Some of the variables listed in Table III. 1 were held constant for the present analysis, but can be 
used if their statistical variation properties are known. As can be seen from the variable definitions 
above, there are motor to motor variations and variations within a single pair of motors. Variations within 
a pair are identified SHFT after the variable name. For example A1SHFT is the within pair variation of 
the burn rate coefficient Al. The analysis consisted of selecting input parameters for 2000 5 -segment 
RSRMV motors or 1000 motor pairs. The variations between motors of a pair were defined for what was 
referred to as the odd and even numbered motors in that pair. Specifically, the odd numbered motor was 
considered the base for the pair and the motor to motor variations were imposed on the even numbered 
motor relative to the odd numbered motor. 
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Table III.l Monte Carlo Analysis Steady State Operation Variables. 


Variable Name 

Variable Definition 

RHO 

propellant Density 

A1 

burning rate coefficient 

A1SHFT 

reference burn rate coefficient shift between motors of a single pair 

N1 

burning rate exponent 

ALPHA 

erosive burning rate coefficient 

BETA 

erosive burning rate coefficient 

ROAL 

oxidizer to aluminum ratio 

DE 

exit plane diameter 

DTI 

initial throat diameter 

THETA 

cant angle of the nozzle 

ALFAN 

nozzle cone half-angle 

XT 

aft end taper dimension 

ZO 

aft end taper dimension 

ZC 

aft end taper dimension 

RONDCN 

nozzle end case out of roundness parameter 

RONDCH 

head end case out of roundness parameter 

RONDGN 

nozzle end grain out of roundness parameter 

RONDGH 

head end grain out of roundness parameter 

EXN 

eccentricity in x at nozzle end 

EYN 

eccentricity in y at nozzle end 

EXH 

eccentricity in x at head end 

EYH 

eccentricity in y at head end 

ALPHAH 

angular orientation of ovality at nozzle end 

ALPHAN 

angular orientation of ovality at head end 

ERREF 

reference throat erosion rate 

TGR 

grain bulk temperature 

TGSHFT 

reference PMBT difference between motors of a single pair 

DO 

outside CP grain diameter 

DI 

initial inside CP grain diameter 

THETAG 

aft end taper angle 

LGCI 

initial length of CP propellant grain 

LGNI 

aft end tapered grain length 

THETCN 

nozzle end closure angle 

THETCH 

head end closure angle 

RC 

propellant star grain outside radius 

FILL 

star valley fillet radius 

RP 

initial radius of truncated star 

RIS 

initial radius at the bottom of the truncated star slot 


It should be noted at this point that it is important have a good baseline model prediction before applying the 
statistical procedures used in the analyses described in this paper. To establish that the nominal values for the 
performance parameters did provide a good baseline using the simplified internal ballistic code described in Refs. 
18 and 19, a design optimization program using the same internal ballistics model was used to select the values for 
the nominal variables which gave the best match to the nominal thrust time trace. This was done using the SLS 
static test data as a reference thrust time trace and selecting the performance parameters to be used with a pattern 
search optimization method. The use of this method to match a given thrust time trace has been documented in 
Refs. 20 and 21. The results of such an analysis are shown in Fig. Ill- 1 . 
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Vacuum Thrust vs. Time 
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Figure III-l. Results from design optimization program to define a baseline motor. 


Based on the Monte Carlo variable selection process, input files for the 2000 individual motors were created. 
These input files were then analyzed sequentially and individually to determine the thrust time trace for each motor. 
The results for the odd and even motors were stored separately and the thrust imbalance versus time was calculated 
for each pair of motors. The thrust imbalance versus time data was analyzed statistically at individual time slices. 
The results of this analysis were the mean and standard deviation for the thrust imbalance at each time slice chosen. 
Based on this data a prediction of the K-sigma limit about a zero mean value envelope was constructed as a function 
time. The same value for K was used for this analysis as was used for the ignition transient imbalance study 
described above. This envelope was then compared to the design thrust imbalance limits for a pair of motors during 
steady state operation. There are two sets of data presented below. One is for the propellant mean bulk temperature 
(PMBT) of 60° F +/- TGSHFT and the other is for a mean value of PMBT selected randomly between 40°F and 90° 
F with the same TGSHFT. The 40° F to 90° F was chosen since it represents the allowable operational range for 
PMBT. The bum rate reference temperature was 60° F for each set of data. Plots for these cases are given Fig. III-2 
shows the Monte Carlo prediction envelope for 1000 motor pairs compared to the design specification envelope. 
The prediction shows values of thmst imbalance which are potentially higher than the design envelope with the peak 
values occurring slightly earlier in time. Fig III-3 sows the same envelopes with the thmst imbalance versus time 
for 1000 motor pairs superimposed. The predicted thmst imbalance versus time curves for the 1000 pairs falls well 
within the both the Monte Carlo prediction and the design envelope. The results for dispersions between motors of 
a single pair about a mean value of 60° F PMBT agree well with the specification limits provided by NASA with 
regard to the magnitude of the thmst imbalance as a function of time. It can be shown that this same agreement, 
indeed the same results can be obtained for dispersions about other values of PMBT, assuming that the bum rate 
reference temperature is the same as the nominal PMBT. The figures also show that the time at which tail off begins 
for both the predicted and design envelopes are also shown to be in good agreement. Fig. III-4 shows the thmst time 
traces for all 2000 individual motors and is an indication of the variation in performance for the overall population 
of motors. Figs. III-5 and III-6 present a comparison of a typical motor pairs thmst imbalance with the envelopes 
and the thmst time traces for each motor in this pair respectively. 
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Thrust Imbalance 
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Figure III-2. Steady state thrust imbalance envelopes 
(PMBT mean = 60° F and K=3.472). 


Vacuum Thrust Imbalance vs. Time for 1000 Motor Pairs 



Figure III-3. +/- Absolute value of the steady state thrust imbalance 
(PMBT mean = 60° F and K=3.472). 
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Thrust Imbalance 


Vacuum Thrust vs. Time for 2000 Motors 



Time 

Figure III-4. Thrust time traces for 2000 motors 
(PMBT mean = 60° F). 


Vacuum Thrust Imbalance vs. Time for Motor Pair 404 



Figure III.5. +/- Absolute value of the thrust imbalance for a single motor pair 
(PMBT mean = 60° F and K=3.472). 
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Vacuum Thrust vs. Time for Motor Pair 404 



Figure III.6. Thrust time traces for a single motor pair 
(PMBT mean = 60° F). 


Figs. III-7 through III- 11 presents the same type of data as Figs. III-2 through III-6, but for the case where the 
PMBT is allowed to vary between 40° and 90° F for the odd motor with the same allowable temperature shift 
between motors of a single pair as was used previously. The reference bum rate temperature is again 60° F. The 
thmst imbalance envelopes are similar during the steady state portion of the bum but differ significantly with respect 
to the tail-off. The time at which the effects of tail-off begin occur are earlier and last for a longer period of time 
than the previous case when the PMBT was 60° F. This is an obvious result, but it indicates that the prediction 
captures accurately the effects of PMBT on the motor bum time. The hotter motors burn out earlier and the colder 
motors burn out later, hence the wider tail-off region and the spreading out in time of the larger thmst imbalances. 
Fig. III-7 shows that on average the web action time occurs earlier in time than is shown in the NASA spec limit 
envelope. And that the maximum thmst imbalance can occur over a wider time range. This is due to the fact that 
the burn rate reference temperature was maintained at 60° F and not changed to the PMBT as was done in 
developing the NASA spec limit. Hence the results shown in Figs. III-7 through III-l 1 accounts for the sensitivity 
of the burn rate to PMBT. 

It should be noted that the overall K- Sigma limit is reduced for this case relative to the one with the temperature 
sensitivity of bum rate removed. This would indicate that the NASA spec limit may be conservative with regard to 
the magnitude of thmst imbalance, but does not account for the difference in web action times for motors operating 
at different values of PMBT which are within the operational limits for the motors. It should be noted that the 
within pair variation of PMBT was the same for each range analyzed. The mean value of the web action time was, 
as would be expected, strongly dependent on the value of PMBT. 
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Thrust Imbalance 


Thrust Imbalance Envelopes 



Figure III-7. Steady state thrust imbalance envelopes 
(PMBT mean = 40° - 90° F and K=3.472). 


Vacuum Thrust Imbalance vs. Time for 1000 Motor Pairs 



Figure III-8. +/- Absolute value of the steady state thrust imbalance 
(PMBT mean = 40° - 90° F and K=3.472). 
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Thrust Imbalance Thrust 
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Figure III-9. Thrust time traces for 2000 motors 
(PMBT mean = 40° - 90° F). 

Vacuum Thrust Imbalance vs. Time for Motor Pair 393 



Figure III-10. +/- Absolute value of the thrust imbalance for a single motor pair 
(PMBT mean = 40° - 90° F and K=3.472). 
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Vacuum Thrust vs. Time for Motor Pair 393 



Figure III-ll. Thrust time traces for a single motor pair 
(PMBT mean = 40° - 90° F). 


To verify that the NASA Spec Limit was independent of PMBT, an analysis using 90° F for mean PMBT with a 
90° F bum rate reference temperature was also analyzed. The results for the thmst imbalance and web action times 
were identical with those above for the 60° F case described above and are shown in Figs. Ill- 12 and III- 13 below. 


Thrust Imbalance Envelopes 



Figure III-12. Steady state thrust imbalance envelopes 
(PMBT mean = 90° F and K=3.472). 
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Vacuum Thrust Imbalance vs. Time for 1000 Motor Pairs 



Figure III-13. +/- Absolute value of the steady state thrust imbalance 
(PMBT mean = 90° F and K=3.472). 


IV. Summary 


Two existing legacy solid rocket motor internal ballistics codes, one for the ignition transient and one 
for the steady state and tail-off periods of operation have been interfaced with a Monte Carlo statistical 
program to predict the thrst imbalance between pairs of solid rocket motors to be used on the SLS vehicle. 
The thrust imbalance envelopes generated by the current analysis are in general agreement with the design 
specification requirements for the SLS vehicle during the ignition transient and the steady state and tail-off 
periods of operation. 
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